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In this paper, we present a theoretical investigation of phot oabsorpt ion spectrum of the newly 
synthesized hydrogenated cluster of aluminium, AI4H6. The calculations are performed within the 
wave-function-based semi-empirical method employing the complete neglect of differential overlap 
(CNDO) model, employing a large-scale configuration interaction (CI) methodology, and our re- 
sults are found to be in very good agreement with the earlier ones obtained from the ab initio 
time-dependent density functional theory (TDDFT). We carefully analyze the many-particle wave 
functions of various excited states up to 8 eV, and find that they are dominated by single-particle 
band to band excitations. This is in sharp contrast to bare aluminium clusters, in general, and AI4, 
in particular, whose optical excitations are plasmonic in nature. We attribute this difference to be 
a consequence of hydrogenation. 

PACS numbers: 36.40.Vz, 31.10.+Z, 31.15.bu, 31.15.vq 



INTRODUCTION 



Because of the future importance of hydrogen as a 
fuel, now-a-days a considerable amount of research ac- 
tivity is taking place in the field of hydrogen storage 
materials p], @|. This is because hydrogen is an extremely 
combustible gas, therefore, the preferred way to store it 
is in physically or chemically adsorbed forms. It is from 
this point of view that a large amount of theoretical and 
experimental research is taking place in the field of metal- 
lic hydrides. Among all possibilities, hydrides of group 
III elements, boron and aluminium, which are called bo- 
ranes and alanes, respectively, are considered to be strong 
candidates for the purpose [3|-[5|. It is with this aim in 
mind, that Li et al.[&\ recently synthesized and studied 
the properties of a new alane, AI4H6. The same group fol- 
lowed up this work, by performing an experimental and 
theoretical study of the electronic structure and prop- 
erties of closo-alanes, AI4H4, Al n H n+ 2, 4 < n < 8@- 
Besides their possible applications, alanes are of inter- 
est from a fundamental point of view as well. Boron 
and aluminium are in the same group of periodic table 
with three valence electrons each, yet their properties are 
very different. In the bulk form boron is a semiconduc- 
tor, while aluminium is a metal. As far as the hydride 
chemistry is concerned, boron exhibits a huge variety of 
boranes with n- vertex polyhedral structures]!, Q , while 
the number and types of alanes is much smaller. The pio- 
neering works of Lipscomb [10] which provide the explana- 
tions of the cage structures of boranes using the concepts 
of three-centered two-electron(3c-2e) bonding, ultimately 
resulted in an ingenious molecular orbital theory known 
as polyhedral skeletal electron pair theory (PSEPT) or 
simply Wade-Mingos rule |lll - [l3| . On the other hand 
only a small number of hydrides of aluminium have 
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been investigated theoretically and experimentally such 
as molecules AffljliHjjj, A1H 3 [HSjAk H 4 jlj Al 2 H 6 [l^ 
[HI, of the general type Al n H n+ 2[25|, |26|, a nd cage-like 
clusters such as Ali3H~ m [27]. Fu et al. [26] in a recent 
comparative study investigated theoretically Al n H n 
clusters, and their borane analogues. Grubisic et al. 
have tried to arrive at a Wade-Mingos type of rule set 
for aluminium hydrides, so that their structures could be 
predicted in a way similar to boranes, though Martinez 
and Alonso[25] have recently raised doubts on this sim- 
ilarity. With this background in mind, the recent study 
of Ai4H6,[6] a material whose structure is consistent with 
Wade-Mingos rule, has again revived the analogy of alu- 
minium hydrides with boranes. In this work we present 
an extensive study of the electronic structure and op- 
tical properties of AI4H6, with the aim of predicting its 
photoabsorption spectrum, and to understand the nature 
of its optical excitations, which can be used for optical 
characterization of this substance. For the purpose, we 
use our recently developed CNDO/INDO Hamiltonian- 
based semiempirical multi-reference singles-doubles con- 
figuration interaction (MRSDCI) methodology described 
elsewhere [28-30]. To benchmark our approach, we also 
apply this approach to A1H molecule and demonstrate 
good agreement with the experimental results. Martinez 
and Alonso[25] calculated the photoabsorption spectra of 
several alanes of the type Al n H n+ 2 including AI4H6, us- 
ing the ab initio time-dependent density functional the- 
ory (TDDFT), and our results are found to be in very 
good agreement with their results [25]. Upon analyzing 
the many-particle wave functions of the excited states 
corresponding to important peaks in the spectrum of 
AI4H6, we conclude that all the way up to 8 eV the 
states correspond to inter-band excitations, and not to 
plasmonic excitations common in bare Al clusters. 

The remainder of this paper is organized as follows. In 
section [III we briefly describe the theoretical methodology 
employed for the present calculations. This is followed by 
the presentation and discussion of our results in section 
IIIII Finally, in section HVl we present our conclusions. 
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II. THEORY 

For our study, we adopted a wave-function based 
electron-correlated approach employing the semi- 
empirical valence-electron CNDO/2 model Hamiltonian 
developed by Pople and coworkers (3ll - [33| . The method- 
ology adopted in this work is discussed in detail in 
our earlier paper s|28l430j. therefore, we present only a 
brief description of it here. As compared to our earlier 
works on boron-based clusters [29], l30l | where we had 
used the INDO model, the choice of CNDO/2 model 
here has to do with the fact that INDO parameters are 
not available for the third row atoms like aluminium 
in the original Pople approach [34]. Our calculations 
are initiated at the Hartree-Fock (HF) level, within 
the CNDO/2 model, using a computer program de- 
veloped recently by us[28]. The CNDO-HF molecular 
orbitals (MOs) thus obtained, are used to transform the 
Hamiltonian from the original atomic-orbital (AO) to 
the MO representation, which is subsequently used in 
the post-HF correlated calculations. The transformed 
CNDO/2 Hamiltonian matrix elements in the MO repre- 
sentation are supplied to the computer program package 
MELD [35], which is used to perform the correlated 
calculations using the multi-reference singles-doubles 
configuration-interaction (MRSDCI) approach. Using 
the ground- and excited-state wave functions obtained 
from the MRSDCI calculations, electric dipole matrix 
elements are computed and subsequently utilized to 
compute the linear absorption spectrum, under the 
electric-dipole approximation, assuming a Lorentzian 
line shape. By analyzing the excited states contributing 
to the peaks of the computed spectrum obtained from 
a given calculation, bigger MRSDCI calculations are 
performed with a larger number of reference states. 
This procedure is repeated until the computed spectrum 
converges within an acceptable tolerance. In the past, 
we have used such an iterative MRSDCI approach on a 
number of conjugated polymers to perform large-scale 
correlated calculations of their linear and nonlinear 
optical spectra[36]. 

III. RESULTS AND DISCUSSION 

In this section we present and discuss our results on 
the electronic structure and optical properties computed 
using our CNDO-MRSDCI approach. However, before 
that, to benchmark our methodology we present and dis- 
cuss our CNDO-CI results on the simplest hydride of 
aluminium, namely A1H. 

A. A1H Molecule 

First we performed the geometry optimization of the 
A1H molecule at the CNDO-HF level, by calculating the 
total energy of the system for different bond lengths, and 



then locating the minimum. Our optimized bond length 
was found to be 1.75 A which is 0.1 A larger than the 
reported experimental value of 1.65 Aflii fig. However, 
our aim here is to compute the optically-active excited 
states of the molecule, and as we will show later that this 
much of difference in the geometry leads to insignificant 
differences in their excitation energies. 



In the CNDO/2 model Al has nine Slater- type basis 
functions (Is, 3p, 5d), while H has only one basis func- 
tion, leading to ten basis functions in all. Being a va- 
lence electron approach, the total number of electrons 
being explicitly considered for A1H within the CNDO/2 
model is four, with Al contributing three electrons, and 
H one. Therefore, with four electrons, and ten basis 
functions full-CI (FCI) calculations for the system are 
feasible, because in the singlet subspace the total num- 
ber of all possible configuration state functions (CSFs) is 
only 825. Thus, the FCI results presented here are exact 
within the chosen model (CNDO/2), and any disagree- 
ments with the experiments will indicate the deficiency 
of the model rather than that of the CI expansion. 



The CNDO-HF calculations predict the Mulliken pop- 
ulation of A1(H) atoms to be +0.24(-0.24), implying 
a partial ionic character to this system, with Al be- 
ing the electron donor, fully consistent with previous 
works [23]. The absorption spectrum of A1H is well-known 
and has been studied extensively over the years both 
experimentally [16] and theoretically fl5l 1171. 1 1811. The 
ground state of the system is classified as while the 

dipole connected excited states are .D 1 S + , 
and E 1 !!. The first excited state A 1 !! corresponds to 
HOMO(#)^LUMO(L) transition which is 5a 2tt in 
nature. The reported experimental value of this transi- 
tion is 2.91 eV, with which our calculations are in good 
agreement with the computed values of 3.11 eV, and 
3.12 eV at the bond lengths of 1.75 A (our optimized 
bond length) and 1.65 A (experimental bond length) 
respectively. This also shows that small differences in 
geometry have negligible impact on the calculated exci- 
tation energy of the lowest excited state. The C and 
D S-bands correspond to Rydberg-type excitations and, 
therefore, cannot be computed using CNDO approach, 
because they require the use of extremely diffused basis 
functions (l5l|. The E 1 !! band corresponds to \H—1 — >• L) 
(4<T — »> 2tt) excitation, and has been measured close to 
6.9 eV[16], while our calculation predicts its value to be 
7.1 eV, thus overestimating it by 0.2 eV. We note that 
quantitatively speaking, the disagreement between our 
results, and the experimental ones, both for the A 1 !!, 
and E 1 !! bands is 0.2 eV, which is quite satisfying, given 
the semiempirical nature of this approach. Having bench- 
marked the CNDO-CI approach for the case of the A1H 
molecule, next we present the results of our calculations 
for A1 4 H 6 . 



B. A1 4 H 6 

We first carried out the geometry optimization of 
AI4H6 within a density functional theory (DFT) based 
approach, using Gaussian03 computational package [37] 
employing B3LYP functional, and 6-311+g(d) basis set. 
The optimized geometry is presented in FigOQ and it is 
virtually identical to that reported by Li et al. [6] As noted 
by Li et al.[6] also, this structure is very similar to the 
computed structure of the borane, B4He[38] which still 
has not been synthesized. We used this DFT optimized 
geometry to perform the calculations of the photoabsorp- 
tion spectrum within our CNDO-CI approach, as well as 
the ab initio TDDFT method. 




Figure 1: Optimized geometry of AI4H6 obtained from 
the ab initio DFT calculations (see text for details). Each 
bond length (in A) has also been indicated. 

Before discussing our calculated spectra, we present 
some of the molecular orbitals (MOs) of the system close 
to the Fermi level. In Fig. [2] MOs obtained from the 
CNDO-HF calculations are presented. Although, here 
we have not presented the ab initio B3LYP MOs which 
were used to perform the TDDFT calculations, we note 
that among the occupied orbitals, H — 3, H — 1, and 
H obtained from the CNDO-HF and the DFT calcula- 
tions were qualitatively very similar. Mulliken charges 
of various atoms indicate Al to be in slightly cationic 
state, while the H atoms carry small negative charges. 
Therefore, we conclude that the bonding in the system is 
largely covalent, with some polar character. The covalent 
nature of the bonding is also obvious from the MO plots, 
which show that there is significant amount of charge 
density in the region between the atoms. 

Next, we present and discuss the linear optical ab- 
sorption spectrum of the system computed using the 
CNDO-CI and the ab initio TDDFT methods. The cal- 
culations were performed using the geometry shown in 
Fig. [TJ which as mentioned above, was obtained at the 
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(a) HOMO-3 (b) HOMO-2 (c) HOMO-1 




(d) HOMO (e) LUMO (f) LUMO+1 (g) LUMO+2 



4 

(h) LUMO+3 

Figure 2: (Color online) Molecular orbitals (iso plots) of 
A1 4 H 6 from HOMO-3 to LUMO+3, obtained from the 
CNDO-HF calculations. 



B3LYP/6-311+g(d) level of theory. For the CNDO-CI 
calculations the total number of valence electrons was 18 
and the total number of basis functions was 42, thus rul- 
ing out the FCI calculations. Therefore, as mentioned 
in section HIl the MRSDCI method was adopted to com- 
pute correlated wave functions and energies of the ground 
and the optically active excited states. The calculations 
were performed without the use of point-group symme- 
try, therefore, the CI expansion became fairly large scale, 
with the largest CI calculation consisting of 1.3 millions 
CSFs. The CI matrix was iteratively diagonalized 35| 
to obtain 60 lowest roots which required several hours of 
CPU time. The dipole moments connecting these excited 
states to the ground state were used to compute the pho- 
toabsorption spectrum presented here. Before discussing 
our final spectrum, in Fig. [3] we present the convergence 
of our optical absorption results with respect to the size 
of the CI matrix, represented in this case by the num- 
ber of reference configurations (N re f) from which singly- 
and doubly- excited CSFs were generated to obtain the 
MRSDCI expansion. For example, in our largest MRS- 
DCI calculation we used N re f=22 which generated 1.3 
million total CSFs. The convergence of our calculations 
with respect to N re f is obvious from Fig. [3j in which 
the spectra corresponding to N re f=14 and N re f=22 are 
virtually indistinguishable. 

In Fig. 2J we present our final linear optical absorp- 
tion spectrum of AI4H6 computed using the CNDO-CI 
method. Before discussing our spectrum in detail, we 
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Figure 3: (Color online) Convergence of the linear 
absorption spectrum of AI4H6 computed using the 
CNDO-MRSDCI method with respect to the number of 
reference states used in MRSDCI calculations. A line 
width of 0.1 eV was used to compute the spectra. 



would like to benchmark it against other calculations to 
ensure its correctness. Therefore, to ascertain the ac- 
curacy of our CNDO-CI calculations, we also performed 
ab initio TDDFT calculation of its lowest photoexcited 
state, at the B3LYP/6-311+g(d) level of theory using the 
Gaussian03 program package[37]. Our CNDO-CI value 
of 2.98 eV (peak I in Figj4j), compares excellently with the 
TDDFT value of 3.03 eV, which gives us confidence as to 
correctness of our calculation. Martinez and Alonso[25j 
used an ab initio TDDFT methodology to compute the 
photoabsorption spectrum of AI4H6 all the way up to 
12 eV (cf Fig. 7 of Ref.flH]). The first peak in the 
TDDFT spectrum of Martinez and Alonso[25] is at an 
energy slightly higher than 3 eV, and their highest peak 
occurs at an energy close to 9 eV, thus making their spec- 
trum slightly blue-shifted compared to ours. However, 
the qualitative nature of their spectrum[25] is quite sim- 
ilar to our spectrum discussed below. 

Next, we examine the nature of excited states corre- 
sponding to the peaks present in our calculated spec- 
trum (cf Fig. The important peaks in the spectrum 
up to the excitation energy «8 eV have been labeled, 
and the many-particle wave functions of these excited 
states, along with their ground-state transition dipole 
moments, are presented in table [U The first peak in 
Fig. [4] occurs at 2.98 eV, and is a relatively weak peak 
corresponding to the across the gap H — >• L transition. 
The next peak (II) located at 3.28 eV is also similarly 
weak and corresponds to H — 1 — >> L excitation. When 
we examine the nature of the many-particle wave func- 
tions (cf Table [I]) of various excited states, we realize 
that: (a) the wave functions of all the states consist pre- 
dominantly of singly-exited CSFs, (b) all the states have 
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Figure 4: Linear optical absorption spectrum of A1 4 H 6 , 
computed using the CNDO-MRSDCI approach. Peaks 
with energies up to 8 eV have been labeled. A line 
width of 0.1 eV was used to compute the spectrum. 



one dominant configuration with magnitude of its coef- 
ficient always greater than 0.8, with some configuration 
mixing with the increasing excitation energy. The high- 
est peak in the spectrum is located at 8.05 eV with the 
dominant CSF being | H — 1 — >• L + 7), with some con- 
tributions from I H L + 9) and | H - 2 L + 3). 
Koutecky and coworkers [39] have formulated a criterion 
according to which if the many particle wave function of 
a given excited state is dominated by one singly-excited 
configuration, it is classified as a normal inter-band ab- 
sorption. On the other hand, if the wave function ex- 
hibits strong mixing of several CSFs with coefficients of 
almost equal magnitude, it is considered to be a plas- 
monic collective excitation [39]. In metallic clusters such 
as AI4H6, to know whether the excited states are plas- 
monic in nature or not, is always of interest. As per the 
Koutecky [39| criterion, the excited states correspond- 
ing to peaks I to VIII do not appear to be of plasmonic 
type. This is to be contrasted with the nature of optical 
absorption in the bare AI4 rhombus cluster which were 
found to be of plasmonic type by Deshpande et al. [40] in 
an ab initio time-dependent local-density approximation 
(TDLDA) calculation. In their work, authors reported 
that low-intensity optical absorption in AI4 begins close 
to 1 eV, and increases in intensity with energy peaking 
around 8 eV. This suggests that in AI4H6, it is the in- 
fluence of hydrogenation which increases the optical gap, 
as well as changes the nature of excited states from plas- 
monic to inter-band excitations. 



IV. CONCLUSIONS 

In conclusion, we have preseneted the linear optical ab- 
sorption spectrum of AI4H6 calculated using a wave func- 
tion based methodology employing the CNDO-CI ap- 
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Table I: Excitation energies and many-particle wave 
functions of excited states corresponding to the peaks in 
the CNDO-CI linear absorption spectrum of AI4H6 (cf. 
Fig. [4j), along with the squares of their dipole coupling 
(/i 2 = J2i \ (f\di\G)\ 2 ) to the ground state. |/) denotes 
the excited state in question, \G), the ground state, and 
di is the z-th Cartesian component of the electric dipole 
operator. In the wave function column, the numbers in 
the parentheses are the CI coefficients of a given 
electronic configuration. Symbols H/L denote 
HOMO/LUMO orbitals. 
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proach, and compared it with that computed using the ab 
initio TDDFT approach. The photoabsorption spectra 
computed using the two approaches are generally in good 
agreement with each other, with the published TDDFT 
spectrum [25] being slightly blue-shifted compared to our 
work. Therefore, it will be of considerable interest to per- 
form optical absorption experiments on AI4H6 cluster, to 
determine the nature of its optical excitations, and shed 
further light on its electronic structure. Our results can 
then be used for optical characterization of this cluster. 
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